Localized mode hybridization by fine tuning of 2D random media 
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We study numerically interaction of spatially localized modes in strongly scattering two-dimensional media. 
We move eigenvalues in the complex plane by changing gradually the index of a single scatterer. When spatial 
and spectral overlap is sufficient, localized states couple and avoided level crossing is observed. We show that 
local manipulation of the disordered structure can couple several localized states to form an extended chain 
of hybridized modes crossing the entire sample, thus changing the nature of certain modes from localized to 
extended in a nominally localized disordered system. We suggest such a chain is the analog in 2D random 
systems of the ID necklace states, the occasional open channels predicted by J.B. Pendry through which the 
light can sneak through an opaque medium. 
© 2013 Optical Society of America 
OCIS codes: 290.4210, 260.2710 



In finite open random media, the properties of wave 
propagation depend on the nature of the modes. When 
scattering is sufficiently strong, modes may be spatially 
localized by disorder [1] if their spatial extension, given 
by the localization length ^, is smaller than the system 
size L. Their linewidth decreases exponentially with dis- 
tance from the sample boundaries [2,3]. In this case, 
wave propagation practically comes to a halt beyond a 
length of a few ^, except at the eigenfrequencies of lo- 
calized states where resonant transmission is maximal in 
a very narrow linewidth. In fact, extended states may 
occur accidentally in the very regime of localization, as 
predicted by J. Pendry [4,5]. These so-called necklace 
states form multipeaked extended states, which span the 
total system. They result from the coupling of several 
modes which overlap both in space and in frequency. 
These modes are rare but they play a dominant role 
in transport since they are spectrally broad. They have 
been observed experimentally by several groups in one- 
dimensional (ID) random systems [6-10]. This picture 
was extended to larger dimensions [4] but the existence 
of such chains of modes in 2D and 3D was never con- 
firmed, probably because of their extreme scarcity. In 
this work, we raise the following question : can a par- 
ticular random structure be designed to form necklace 
states instead of looking for such rare events in a collec- 
tion of random realizations ? It was shown that intro- 
ducing correlations in disordered structures may force 
the modes to extend, destroying at the same time local- 
ization [11,12] or conversely, improve localization [13,14]. 
However introducing correlation is a strong perturbation 
which affects fundamentally the transport properties of 
the system. 

In this Letter, we show that manipulation of the ran- 
dom structure at the scale of a single scatterer is suffi- 
cient to change drastically the nature of the modes of 
a strongly disordered system. We numerically study the 



spectral and spatial evolutions of localized modes in fi- 
nite 2D random systems as the refractive index of a sin- 
gle scatterer is progressively changed. We observe mode 
crossing and anticrossing as modes overlap spectrally 
and spatially. These modes hybridize to form double- 
peak modes. We demonstrate that by changing the scat- 
tering index of only 2 scatterers, it is possible to couple 
several localized modes and tailor an extended mode, 
much like ID necklace states, without changing the na- 
ture of the random system. Transmission measurements 
show that transport is enhanced through such an open 
channel. 

The system is an open 2D finite structure made of 896 
circular dielectric scatterers with high refractive index 
{ric = 2) embedded in a host material of index unity 
(vacuum) [15]. We work in the range [434. 1,441. 6] nm 
([6.8, 6.9]10^^IIz) with scatterers radius r = 60 nm 
and surface filling fraction <l> = 40%. The system area 
is L X L = bfim x 5/im. Estimations of the localiza- 
tion length gives ^ = Ifim <C L. We solve Maxwell's 
equations using a finite element method [16] to ob- 
tain the eigenstates and their complex eigenfrequencies, 
u = -\- \v" in the TM configuration (electric field E 
perpendicular to the 2D propagation plane). All modes 
in this frequency range are found spatially localized in 
agreement with the estimated value of ^. 

The refractive index of a single scatterer is varied 
from Us = 1.05, to Uq = 3.2, by steps dug = 0.05. For 
each value of , all complex eigenfrequencies and eigen- 
functions in the frequency range considered are com- 
puted. As they are localized away from the scatterer, 
most of the modes are insensitive to this local pertur- 
bation. A few modes however are shifted in frequency. 
The displacement of the real and imaginary parts of the 
eigenfrequencies can be followed " adiabatically" as the 
refractive index increases. This is seen in Fig. 1 where 



a subset of three eigenvalues, Vj = v'. 
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are plotted for each value of (named Mode^^fi 2,3] for 
all values of ng). 
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Fig. 2. (Color online) Spatial distribution of the field am- 
plitude for Model and Mode2 of Fig.l for three different 
values of the refractive index of the scatterer located at 
the white spot. Media 1 displays the evolution of the 
mode profiles as Ug increases from 1.05 to 2.40. 



Fig. 1. (Color online) Monitoring of (a) the real part 
and (b) the imaginary part of three localized states as 
the refractive index, n^, of a particular scatterer (shown 
in Fig. 2) is varied from 1.05 to 2.70. Y-axis has been 
split in (b) to span a broad interval of 

When Us increases, the real frequency of Modei is 
seen to decrease. As Modei and Mode2 get close spec- 
trally, level repulsion is observed for the real part of the 
frequency, u'. Simultaneously, the imaginary parts, 
cross. 

The corresponding spatial distributions of their eigen- 
functions are shown in Medial and Fig. 2 as Us increases 
from 1.05 to 2.40. As Modei and Mode2 get closer to the 
point of avoided crossing at Ug = 2.25, mode hybridiza- 
tion occurs. Two double-peaked modes are formed at 439 
nm and 441 nm with a spectral width of 0.85nm. These 
two modes have identical amplitude maps but different 
phase maps (not shown). Beyond the avoided crossing 
point, the spatial patterns are exchanged. We note that 
the spatial distribution of Modei remains unaffected as 
long as it is spectrally away from other modes. This indi- 
cates that the modal deformation is mainly due to mode 
interaction rather than to the small local perturbation. 

Modei and Mode2 are modes well localized inside 
the system, which have small leakage through the open 



boundaries. Hence, like in purely Hamiltonian systems 
[17], mode interaction leads to level repulsion (anticross- 
ing). If the system is not conservative, the energies be- 
come complex. In this case, depending on the interac- 
tion, crossing as well as anticrossing of energies can oc- 
cur [18]. This is true for quantum systems as well as 
for classical waves [19]. In open random media, leakage 
of modes localized close to the open boundaries can be- 
come significant. Hence, crossing can also be expected. 
This is illustrated by the interaction between Modei and 
Modes which is a mode close to a boundary (not shown). 
Their real parts, ^ cross while repulsion is seen for the 
imaginary part, (see Fig. 1). 

The double-peak hybridized modes seen in Fig. 2 for 
Us = 2.25 present an unexpected spatial profile and ex- 
tension in the regime of Anderson localization, where 
nominal modes are exponentially localized single-peak. 
It was therefore tempting to manipulate further the lo- 
calized modes to increase the number of hybridized lo- 
calized states and try form a chain that connects one 
end of the system to the other. Such an extended hy- 
bridized mode would provide an open channel through 
which light can traverse this nominally localized random 
system, changing radically its transport properties. We 
identified modes localized on both sides of the hybridized 
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mode of Fig. 2 (n5=2.25), which could open completely 
the channel from left to right. We found two scatterers 
which couple all four modes. It is remarkable that such 
an extended mode may coexist with localized modes. It 
is important to realize that the perturbation is small and 
that most of the modes are virtually unchanged, in con- 
trast to e.g. correlated disorder. 

We probe the predominance of this beaded mode in 
transport by comparing in Fig. 3 transmission spec- 
tra through the perturbed and unperturbed systems. A 
plane wave is incident on the left side of the system and 
the total outgoing intensity at the right side is recorded 
as a function of the wavelength of the incident wave. 
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Fig. 3. (Color online) (a) Hybridization of 4 localized 
states by setting the index of refraction to 2.9 and 2.8 for 
left and right scatterers shown as white spots, all other 
scatters set at Uc = 2.0. (b) Transmittance spectra for 
plane wave impinging from left on perturbed (full line) 
and unperturbed (dotted hue) systems, (c) Field am- 
plitude spatial distribution for wave impinging on per- 
turbed system at A = 439.60nm. (d) Field amplitude 
spatial distribution for wave impinging on unperturbed 
system at A = 440.63nm. 

The two narrow peaks of the unperturbed system at 
438 nm and 441 nm correspond to the resonant excita- 
tion of Mode 1 and Mode 2 respectively. In contrast, the 
perturbed system shows a higher and broader transmis- 
sion near 439 nm corresponding to the excitation of the 
hybridized mode as shown in Fig. 3. Hence, a channel 
has been open with enhanced broad band transmission. 
Further investigation is necessary to confirm the predic- 
tion that such necklace states dominate the conductance 
even in 2D random systems [21]. 

In conclusion, level repulsion was observed between 
two localized modes of a disordered system as the index 
of refraction of a single scatterer was increased. For these 
two modes, hybridization is observed at the anticrossing 
point with the formation of symmetric and antisymmet- 



ric double peaked modes. Level crossing can also occur 
when one mode is strongly leaking out of the sample. 
We show that the random structure can be manipulated 
further to create a chain of hybridized localized modes, 
extended from one end of the sample to the other. The 
possibility of forming such open channels is of extreme 
interest. First it confirms that there is an analog in 2D 
random systems of the necklace states discussed in ID 
random media. In ID systems, J. Pendry predicted that 
these extended modes dominate the statistics of conduc- 
tance. This prediction is not obvious for larger dimen- 
sions and needs more investigation. It is tempting how- 
ever to make a connection between open transmission 
channels [20] and necklace structures as this was sug- 
gested by Pendry in the ID case [21], a question to be 
addressed in future work. 
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